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Surveys with the Spitzer and Herschel space observatories are now enabling the discovery and characterization of large 
samples of protostars in nearby molecular clouds, providing the observational basis for a detailed understanding of star 
formation in diverse environments. We are pursuing this goal with the Herschel Orion Protostar Survey (HOPS), which 
targets 328 Spitzer-identified protostars in the Orion molecular clouds, the largest star-forming region in the nearest 500 pc. 
The sample encompasses all phases of protostellar evolution and a wide range of formation environments, from dense 
clusters to relative isolation. With a grid of radiative transfer models, we fit the 1-870 /im spectral energy distributions 
xj^ , (SEDs) of the protostars to estimate their envelope densities, cavity opening angles, inclinations, and total luminosities. 

After correcting the bolometric luminosities and temperatures of the sources for foreground extinction and inclination, we 
find a spread of several orders of magnitude in luminosity at all evolutionary states, a constant median luminosity over the 
more evolved stages, and a possible deficit of high-inclination, rapidly infalling envelopes among the Spitzer-identified 
sample. We have detected over 100 new sources in the Herschel images; some of them may fill this deficit. We also 
report results from modeling the pre- and post-outburst 1-870 /im SEDs of V2775 Ori (HOPS 223), a known FU Orionis 
outburster in the sample. It is the least luminous FU Ori star with a protostellar envelope. 
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t— I \ 1 Introduction pleted 200 hour Herschel Open Time Key Program to ob- 
■ serve most of the Orion protostars with the PACS instru- 
^7"! \ Nearly a decade after the launch of the Spitzer Space Tele- ment on board Herschel. This includes imaging them at 70 
t> ■ scope, thousands of young stellar objects (YSOs) have been and 160 \im and obtaining spectroscopy of a subset of them 
identified in nearby star-forming regions. YSOs can be clas- from 55 to 200 ^m. Supplementing the Herschel program, 
sifted on the basis of their Spitzer mid-infrared colors as we have Hubble Space Telescope and ground-based near- 
er \ protostars, which still retain their nascent envelopes, or as IR imaging, ground-based near-IR spectra, existing Spitzer 
young stars with only disks (Allen et al. 12004b . With the imaging and spectra, ground-based sub-mm line and contin- 
Herschel Space Observatory, Spitzer-identified protostars uum imaging, and photometry from 2MASS and WISeQ 
can be observed in the far infrared, where their spectral en- Here we describe the SED modeling of the HOPS sam- 
ergy distributions (SEDs) peak, allowing an improved es- pie, present an inclination-corrected BLT (bolometric lumi- 
timation of their luminosities and envelope densities. The nosity and temperature) diagram for Orion, and discuss an 
distributions of these quantities for large populations of pro- FU Ori outburster in the sample, 
tostars in diverse environments are a key tool in constraining 
theories of star formation (Dunham et al. |2010t . 

At 420 pc, the Orion molecular clouds are home to 3480 
Spitzer-identified YSOs, 488 of which are likely protostars 

(Megeath et al. [2013, the most of any star-forming region The HOPS P ro 8 ram ima 8 ed 328 of the Spitzer-identified 

in the nearest 500 pc. To better characterize the Orion pro- ° non Protostars using PACS at 70 and 160 /im in 108 dis- 
tostars, we are conducting HOPS, the Herschel Orion Pro- i The Two Micron A11 sky Survey (2M ASS) and the Wide-field Infrared 

tOStar Survey. The cornerstone of HOPS IS a recently Com- Survey Explorer (WISE) generated point-source catalogs for the entire sky 

at near- and mid-infrared wavelengths, respectively, that can be accessed at 



2 Observations 
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tinct square fields of 5' or 8' on a side. Figure Q] shows the 
locations of the observed protostars. The FWHM of point 
sources are 5.2" at 70 /mi and 12" at 160 /tin. 

For photometry, we measure the flux densities in small 
apertures designed to reduce the influence of the bright and 
non-uniform extended emission in Orion. Our apertures are 
9.6" in radius at 70 yum and 12.8" in radius at 160 /iin, with 
subtraction of the median signal in a background annulus 
extending from the aperture limit to twice that value in both 
channels. The results were divided by measurements of the 
encircled energy fractions in these apertures, with an addi- 
tional correction for the fact that our close-in sky subtrac- 
tion removes 3-4% of the flux density in each point-spread 
function. These aperture corrections are 0.733 at 70 /im and 
0.660 at 160 /im. The calibration accuracy of the PACS pho- 
tometry is 5%. The 70 /im flux densities range from 8 mJy to 
930 Jy, while the 160 /im flux densities range from 1 1 mJy 
to 950 Jy. The full table of HOPS 70 and 160 /mi photome- 
try will appear in a subsequent publication. 

The remaining photometry for the SEDs comes from 
2MASS, Spitzer IRAC and MIPS photometry (Megeath et 
al. 12012b . our Spitzer/IRS spectroscopy campaign (C. A. 
Poteet et al. 2013, in prep.), and our APEX sub-mm imaging 
of Orion (T. Stanke et al. 2013, in prep.). We also include 
100 /im PACS photometry from the Herschel public archive. 
For SED fitting, the IRS spectra are converted into discrete 
flux densities at 9.7 /im to incorporate the depth of the sili- 
cate feature there, and at 12, 18, and 32 /im to constrain the 
mid-IR continuum. 

3 SED Fitting 

To fit the SEDs, we constructed a grid of radiative transfer 
models generated with the Monte Carlo code of Whitney 
et al. (2003). The code models protostars as a central lumi- 
nosity source, a luminous accretion disk with a power-law 
radial density profile and scale height, an envelope with a 
density that follows the rotating collapse solution of Tere- 
bey, Shu, & Cassen (119841) . and a bipolar envelope cavity. 
The SEDs of protostars are most sensitive to the total lu- 
minosity (stellar plus accretion), the envelope density, and 
the cavity opening angle. Reprocessing of the stellar and ac- 
cretion photons by the envelope washes out any distinction 
between stellar and accretion photons; hence, we use one 
central source with mass 0.5 M Q , radius 2.09 Rq, and ef- 
fective temperature 4000 K, varying the disk and envelope 
properties to create the grid. We adopt a distance of 420 pc. 

The HOPS grid was first described in Ali et al. d20T7jb . 
The current version consists of 3600 models in which pa- 
rameters of interest in fitting protostellar envelopes are var- 
ied. There are twenty envelope infall rates that range from 
to 10 -3 M Q yr _1 , five cavity opening angles that range 
from 5° to 45°, three disk radii of 5, 100, and 500 AU, three 
envelope radii of 5000, 10,000, and 15,000 AU, and four 
disk accretion rates ranging from 10~ 8 to 10~ 5 M yr _1 . 
The disk accretion rates effectively set the luminosity of the 
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Fig. 1 Extinction map of the Orion molecular clouds; circles 
mark the Spitzer-identified protostars observed by PACS. The 
HOPS sample was required to have 24 /im photometry to ensure 
a more reliable identification of protostars. Since the center of the 
Orion Nebula was saturated in the Spitzer 24 /im survey of Orion, 
the HOPS sample does not include sources in its brightest regions. 



system, where the central star always contributes 1 Lq and 
accretion contributes the rest. The lowest accretion rate cor- 
responds to a total luminosity of 1.1 Lq, and the highest 
accretion rate corresponds to 56.3 Lq. Each model can be 
viewed from one of ten inclination angles, resulting in a grid 
of 36,000 SEDs. Each SED can be altered in two ways to 
improve the fit. It can be multiplied by a constant, which al- 
lows the best-fit luminosity to fall along a continuum instead 
of four discrete values, and it can be reddened by foreground 
dust up to Ak = 3. A future version of the grid will include 
less luminous central stars, reducing the dependence on lu- 
minosity scaling to model low-luminosity objects. 

4 Bolometric Luminosities and 
Temperatures: Inclination Dependence 

For each protostar, we calculate the bolometric luminosity 
and temperature (BET properties; Lboi and Tboi) of its best- 
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fit model. Myers & Ladd ( 119931 ) defined T bo i to be the tem- 
perature of a blackbody with the same mean frequency as 
the SED. The bolometric temperature increases as the enve- 
lope evolves; Chen et al. dl995b found that T bo i < 70 K for 
the Class sources (deeply embedded protostars), 70 K < 
Tboi < 650 K for the Class I sources (protostars), and 
650 K < Tboi < 2800 K for the Class II sources (young 
stars with disks). To examine changes in protostellar prop- 
erties within Class I, we further divide Class I into early and 
late stages, with the transition at Tboi = 210 K. 

We modify the original definition of Tboi to include only 
the thermal component of the best-fit model SED in the cal- 
culation. (The Whitney et al. code reports the independent 
SEDs of direct light, scattered light, and thermal emission, 
which together sum to the total.) Including only the ther- 
mal emission is more representative of the envelope density; 
protostars can have substantial scattered light in their SEDs 
that shifts Tboi to misleadingly large values. 

Measuring the best-fit model allows corrections to be 
made for incomplete wavelength coverage and the effects 
of foreground reddening and inclination. Foreground red- 
dening shifts the SEDs to lower Tboi and Tboi- An edge-on 
source will yield a lower Tboi and Tboi than the source in- 
trinsically has, while a pole-on source will yield systemat- 
ically high values for both. To correct for foreground red- 
dening, we simply deredden the SED by the best-fit Ak- To 
correct for inclination, we first find the inclination-averaged 
SED of the best-fit model by integrating its flux density over 
all inclination angles and then measure the BET properties 
of this SED, denoted (Tboi) and (Tboi). 

Figure [2] shows the bolometric luminosities and tem- 
peratures of the HOPS sources with correction for fore- 
ground reddening (top) and with an additional correction 
for inclination (bottom). While 26% of the protostars have 
Tboi < 70 K, only 3% have (T bo i) < 70 K, implying 
that sources labeled "Class 0" based on their Tboi alone 
may be more evolved protostars viewed at high inclinations, 
through the denser regions of the envelope or the disk. Ta- 
bles Q] and |2] show statistics for the usual bolometric tem- 
perature categories before and after the inclination correc- 
tion. In both cases, we see a large spread in luminosities 
at each bolometric temperature and a decreasing envelope 
infall rate with increasing bolometric temperature, although 
the latter effect is more pronounced after the inclination cor- 
rection. After the correction, there is a constant median lu- 
minosity for (Tboi) > 70 K. The average inclination an- 
gle for (Tboi) > 70 K is ~ 60°, as expected for a ran- 
dom distribution of orientations. In contrast, sources with 
(Tboi) < 70 K have systematically lower inclinations be- 
cause many of the sources with Tboi < 70 K are fit to 
higher-inclination models with (Tboi) > 70 K. Given the 
strong inclination dependence of Tboi, we are investigating 
alternative diagnostics of envelope evolution such as (Tboi)- 

Rapidly infalling envelopes at high inclinations seem to 
be missing from the sample of Spitzer-identified Orion pro- 
tostars. However, in our PACS 70 /im images, we detect 130 
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Fig. 2 BLT diagrams for 328 Spitzer-identified protostars in 
Orion. Top: Extinction-corrected bolometric luminosity versus 
extinction-corrected bolometric temperature. Dashed lines divide 
the sample into, from left to right, Class II (crosses), late Class I 
(triangles), early Class I (squares), and Class (circles). Bottom: 
The same quantities corrected for inclination. Each source retains 
the same symbol as above, indicating how the sources are reclas- 
sified when the dependence on inclination is removed. 

Table 1 Protostellar statistics before inclination correction. The 
last three columns report medians for each subsample. 

T bo l L Menv i 

(K) Number (L ) (M yr" 1 ) (°) 

0-70 85 TO 7.5 x 10 _b ST~ 

70-210 109 2.0 2.5 x 10" 6 63 

210-650 127 2.6 2.5 x 10~ 7 57 



Table 2 Protostellar statistics after inclination correction. The 
last three columns report medians for each subsample. 

(Tboi) L M cnv i 

(K) Number (L ) (M yr" 1 ) Q 

0-70 10 L2 1.0 x 10" 4 25~ 

70-210 123 2.6 7.5 x 10~ 6 63 

210-650 186 2.6 5.0 x 10~ 7 63 
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point sources that were not classified as protostars by the 
Spitzer study of Megeath et al. (120121 1 because at 24 /im they 
were either undetected or fainter than the selection criterion. 
While many of these are likely extragalactic and other con- 
taminants, a fraction of them have colors consistent with 
edge-on, high-infall-rate protostars, potentially contributing 
to the (Tboi) < 70 K region of the BLT diagram. These 
objects, the PACS Bright Red Sources (PBRS), will be dis- 
cussed by A. M. Stutz et al. in a forthcoming publication. 

5 A Low-Luminosity FU Orionis Object 

Each forming star may experience several FU Orionis out- 
bursts, which are up to century-long increases in luminos- 
ity by up to three orders of magnitude due to enhanced ac- 
cretion from the circumstellar disk (Hartmann & Kenyon 
119961) . These events may be an important factor in the broad 
luminosity distributions of protostars (Dunham et al. 1201 Oi l. 
In their near-infrared monitoring of L 1641, Caratti o Garatti 
et al. ( 120111 1 discovered that HOPS 223 ([CTF93] 216-2; 
V2775 Ori) was several magnitudes brighter in 2010 than 
when it had been observed by 2MASS in 1998. Fischer et 
al. d20 1 2b presented a comprehensive study of its variable 
SED and multiwavelength images. Spitzer data constrain 
the main rise in luminosity to have occurred between 2005 
April 2 and 2007 March 12, over which time the source be- 
came 2.8 mag brighter at 24 /im. By 2008 November 27, the 
source had dimmed by 45% at 24 /xm. Since then, A'-band 
imaging has showed an apparent decline by 0.25 mag, al- 
though further monitoring is required to determine whether 
this is permanent. 

While not all outbursting protostars are members of the 
FU Ori class, near-infrared spectra of HOPS 223 obtained 
by Caratti o Garatti et al. in 2009 and 2010 and by us in 
2011 and 2012 have the distinctive appearance of FU Ori 
sources: no emission lines, broad H2O absorption, and deep 
absorption in the CO 2-0 rovibrational features. The source 
also has deep blueshifted absorption at He I A10830. This 
likely probes a disk wind (Connelley & Greene [20l0jl and 
is the strongest tracer of outflow activity in near-infrared 
moderate-resolution spectra of FU Ori sources. 

HOPS 223 is among the least luminous FU Ori sources, 
with extinction-corrected luminosities of ~ 51 Lq at its 
2007 maximum and 28 Lq when the FU Ori-like near-IR 
spectra were subsequently acquired. This makes it less lu- 
minous than any FU Ori source with a documented outburst 
except for V2493 Cyg (Miller et al. 120111 ). which has a lu- 
minosity of <~ 12 Lq. HOPS 223, however, has a detectable 
envelope, with 0.09 M Q inside of 10,000 AU according to 
our modeling, while V2493 Cyg does not. These objects 
are evidence that the FU Ori phenomenon is important over 
a wide range of evolutionary states and is responsible for 
some of the scatter in BLT diagrams. The estimated accre- 
tion rate of 0.9 - 1.9 x 10~ 5 M yr _1 is toward the low 
end of the range expected for FU Ori stars. As the source 
fades, it may soon develop an emission-line spectrum typi- 



cal of normal accreting young stars (Baraffe, Vorobyov, & 
ChabrierEoHll. 



6 Conclusions 

As part of the Herschel Orion Protostar Survey, we have 
constructed the 1-870 /.tin SEDs of 328 Spitzer-identified 
protostars in the Orion molecular clouds and automatically 
fit them with a grid of radiative transfer models. We plot the 
bolometric luminosities and temperatures of the protostars, 
using only the thermal emission to calculate the bolomet- 
ric temperature. After correcting for foreground extinction 
and inclination, we find a spread in luminosities of several 
orders of magnitude at all bolometric temperatures, a con- 
stant median luminosity over the more evolved stages, and 
a possible deficit of high-inclination, rapidly infalling en- 
velopes in the Spitzer-selected sample. The extremely red 
point sources detected by Herschel but not Spitzer may be 
these missing young protostars. One of the HOPS targets, 
HOPS 223, is a low-luminosity FU Ori outburster. 
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